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ABSTRACT 

We investigate the importance of the shielding of chemical photorates by molecu- 
lar hydrogen photodissociation lines and the carbon photoionization continuum deep 
within models of photon dominated regions. In particular, the photodissociation of 
N2 and CN are significantly shielded by the H2 photodissociation line spectrum. We 
model this by switching off the photodissociation channels for these species behind the 
HHH2 transition. We also model the shielding effect of the carbon photoionization 
continuum as an attenuation of the incident radiation field shortwards of 1102 A. Us- 
ing recent line and continuum cross section data, we present calculations of the direct 
and cosmic ray induced photorates for a range of species, as well as optically thick 
shielding factors for the carbon continuum. Applying these to a time dependent PDR 
model we see enrichments in the abundances of N2, N2H + , NH3 and CN by factors of 
~ 3 — 100 in the extinction band A„ = 2.0 to A^ = 4.0 for a range of environments. 
While the precise quantitative results of this study are limited by the simplicity of 
our model, they highlight the importance of these mutual shielding effects, neither of 
which has been discussed in recent models. 

Key words: astrochemistry - ISM: molecules 



1 INTRODUCTION 

Photodissociation or Photon Dominated Regions (PDRs) 
are pervasive, ubiquitous components of the interstellar 
medium, being the interface between the partially ion- 
ized/atomic component and the neutral, dense, molecular 
component. Examples of PDRs include the boundaries of 
molecular clouds subjected to the interstellar radiation field, 
protoplanetary discs irradiated by accreting protostar and 
the interstellar medium of starburst galaxies illuminated 
by clusters of massive star formation. These interfaces are 
chemically rich and they often dominate the molecular emis- 
sion from the host regions, due to the relatively high exci- 
tation of the molecular gas. Moreover, as earlier, simplistic 
assumptions about the geometries of these regions gives way 
to more realistic, complex morphologies, it would seem that 
a larger proportion of interstellar gas is contained within 
PDRs. 

However, the physical and chemical structure of PDRs 
is complex and only partially described by existing mod- 
els. The thermal balance in PDRs is determined by heat- 
ing by ambient and stellar far-ultraviolet radiation (6 < 
hv < 13.6eV) and cooling by atomi c /molecular line emis- 
sion and continuum emission by dust l|Hollenbach fc Tielend 
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ll999tlSternberell2004h . This in turn depends on the chem- 
istry. The transition from atomic to molecular gas is defined 
by the radiative transfer, coupled to the photochemistry, 
which entails complex self- and mutual-shielding effects. 

It is usually reasonable to assume that H2, CO and its 
isotopologues are the only molecular species that have large 
enough abundances to result in significant absorption of the 
incident radiation field. The primary consequence of this is 
the increasing self-shielding of the photoreaction as the tran- 
sition becomes more optically thick. In addition, where there 
is overlap between the absorption lines of different species, 
mutual shielding will occur. Also, one must consider the con- 
tributions to the UV opacity from the carbon ionization con- 
tinuum (for A < 1102 A) and absorption by dust. This leads 
to a range of possible mutual and self-shielding effects as 
listed below: 

(i) H2 self-shielding, 

(ii) CO self-shielding, 

(iii) H2, CO and CO isotopologue mutual-shielding, 

(iv) Mutual shielding by coincident Lyman transition 
lines of Hi and Lyman and Werner transition lines of H2 , 

(v) Shielding of all species by dust absorption, 

(vi) Shielding of H2 and CO by the carbon ionization con- 
tinuum, and 

(vii) Shielding of other species by all of the above pro- 
cesses. 
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Of these, the cross section for attenuation by the dust 
continuum is approximately constant over the wavelength 
ranges of individual photoreaction channels and so is sepa- 
rable from the other shielding terms. In practice, the theo- 
retical studies are somewhat limited by the paucity of labo- 
ratory and theoretical data for molecules other than H2 and 
CO 

iTielens fc Hollenbachl (|l985h developed a PDR model to 
study the chemistry and thermal balance of gas with den- 
sity 10 3 cm -3 < n < 10 6 cm -3 illuminated by far-ultraviolet 
fluxes x = 10 3 ~ 10 6 times more intense than the inter- 
stellar radiation field. They identified the presence of sev- 
eral physically distinct regions within PDRs, including an 
Hi— >ti2 transition zone at A„ < 1 and a complex, warm 
C + — s-C— s-CO transition zone at A„ ~ 3 — 5. Their model, 
which included carbon continuum self-shielding, yielded Cl 
column densities of 4.5 — 16.0 x 10 18 cm -2 at w hich levels 
the co ntinuum is very optically thick. Similarly, iRawlingi 
( 1988) developed a photochemical model of the (dust-free) 
PDR in the vicinity of a nova, which concentrated on the 
mutual shielding between H2, CO and Cl. 

Subsequent studies, in cluding the PDR compari- 
son/benchmarking study of iRollig et al.l ((2007) , although 
more sophisticated than the earlier models, and employing 
accurate photoreaction data for CO and H2, essentially con- 
firmed the basic structure with the C + — >C— >-CO transition 
occurring at A„ ~ 2 — 5 for n~ 10 3 cm -3 and x ~ 10. These 
various models take careful account of the overlap between 
the CO and H2 absorption lines and include the dust con- 
tinuum, but do not make specific reference to shielding by 
the Cl ionization continuum. 

With the availability of reliable wavelength and os- 
cillator strength data for the photodissociation of CO, 
a much more detailed and accurate treatment of the 
self/mutual shielding of H2 and CO has been possible. 
Ivan Dishoeck fc Blac k ( 1988) considered the photodissoci- 
ation of CO and its main isotopologues together with H and 
H2, including self- and mutual-shielding effects. They also 
included shielding by dust and referred to the attenuation 
by the carbon continuum. As shown in that paper (e.g. see 
Fig. 4), when the CO, Hi and H2 absorption lines become 
very optically thick, they can effectively blanket a signifi- 
cant fraction of the continuum in the 912 A < A < 1079 A 
wavelength range, with implications for the photodissocia- 
tion ratesj)f_other speci es . 

IVisser et all (1200911 included self-shielding, mutual- 
shielding and shielding by H and H2 to calculate the pho- 
todissociation rates of CO and its isotopologues, including 
rarer variants. For these species only, the photodissociation 
rate coefficient can be written as: 

k = x k O [N(H 2 ), N(CO), N(H), 4>] exp(- 7 A,), (1) 

where ko is the unattenuated rate, O is the shielding 
function, exp(— *yA v ) is the (separable) dust extinction term, 
and (p represents other factors, such as the Doppler widths 
of the absorption lines, the excitation temperatures for CO 
and H2, and ratios of the column densities of the CO iso- 
topologues. As with previous studies the line and continuum 
attenuation component s were treate d as s eparable. How- 
ever, as pointed out by Ivan Dishoeckl |l988l) the Ai, depen- 
dence does not include the attenuation at A < 1102 A by 
H2 line/continuum and Cl continuum absorption which will 



be particularly important for species like CN and N2, which 
photodissociate at the shortest wavelengths. 

Thus, the earliest studies of PDRs (e.g. 
ITielens fc Hollenbachl [l985) - although lacking the de- 
tails of more recent models - attempted to deal with 
most of the shielding effects listed above. In subsequent 
studies all discussion of the role of the carbon ionization 
continuum is seemingly absent. This is presumably on the 
assumption - which may or may not be valid - that dust 
extinction dominates. However, the cross section for carbon 
photoionization is 1.6 x 10~ 17 cm 2 . It is fairly reasonable to 
assume that column densities of atomic carbon in excess 
of 10 17 cm~ 2 will pe rtain deep within PDRs (e.g. Fig 11 
of |V isser et al. 2009) and so this opacity source must be 
considered. 

In addition, the more recent PDR models have tended 
to concentrate on the carbon budget, while the effects of 
mutual shielding by H2 and CO on other molecules have 
been largely overlooked. In this study we consider two of 
these effects; the mutual shielding by: 

(i) molecular hydrogen lines, and 

(ii) the carbon ionization continuum. 

Critically, there is extensive overlap of both of these with 
the photodissociation cross sections of CN and N2. The mu- 
tual shielding effects therefore have the potential to pro- 
mote nitrogen-based chemistry at extinctions much shal- 
lower than currently seen in simulations. This paper presents 
shielding factors for a range of species due to these effects 
and investigates their influence on the chemical structure. 

Our calculation of the degree of mutual shielding by H2 
dissociation and carbon ionization are presented in Section 
[2] These shielding factors are then implemented in a simple 
PDR model, described in Section [3] with results presented 
in Section [3] We close with a discussion of our results and 
concluding remarks in Section [5] 



2 PHOTOREACTION RATES 

In PDRs, the shielding of the incident radiation field by 
molecules is typically a combination of line and continuum 
processes leading to a complex radiative transfer problem. 
Most models include self shielding by CO and H2 by adding 
multiplicative shielding factors to their photorates as de- 
scribed above. We extend this procedure to investigate mu- 
tual shielding by the H2 line and Cl ionization continuum 
absorption. In this section we detail our method for cal- 
culating photoreaction rates and shielding factors for these 
species. 

2.1 Direct Photorates 

It is usually assumed that the line and continuum contribu - 
tions to the photorates are separable (|van Dishoeckl Fl988). 
Unshielded rate coefficients, k p d, can be calculated using the 
follo wing equations for line and continuum photoprocesses 
(e.g. Ivan Dishoeckl Il988h : 

k c p T= / o a(A)/(A)dAs-\ (2) 
J912A 
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Table 1. Photoioniz ation Rates , Carb on Continuum Shielding Factors, Cosmic Ray Induced Pho toreactio n Efficiencies an d Cross Section 
Data References; Ivan DishoedJ Il988l ( http : //home.strw.leidenu niv.nl/~ewine/photo/ vD88). iBTanscombl [l963 (B62), iBarsuhnl Il977l 
(B77l. lBlack fc Dalgarnolll977l (BD77) and lRawIings et al.|[l99i (RDB93). 



Species Photorate / s 



Maximal Carbon 
Shielding Factor, Si 



CR Induced Photoreaction 
Efficiency, pjy,/ 



Reference 



c 


3.15 x 10~ 10 


0.000 


264.8 


vD88 


C" 


3.31 x icr 08 


0.994 


1037.8 


B62 


CH 


7.63 x 10- 10 


0.142 


547.7 


vD88 


CH 2 


3.58 x 1CT 10 


0.449 


300.0 


B77 


CH 3 


4.93 x lO" 10 


0.600 


509.5 


BD77 


NH 


1.50 x icr 11 


0.000 


22.0 


B77 


NH 2 


1.73 x icr 10 


0.000 


149.3 


RDB93 


NH 3 


2.82 x lO" 10 


0.298 


230.5 


vD88 


O2 


7.65 x lCT 11 


0.000 


57.7 


vD88 


OH 


2.00 x 10" 11 


0.000 


26.6 


B77 


H 2 


3.44 x lO" 11 


0.000 


24.5 


vD88 


Na 


1.54 x lO" 11 


0.829 


9.1 


vD88 


HCO 


4.81 x ID" 10 


0.590 


499.2 


RDB93 


H 2 CO 


4.82 x 10" 10 


0.215 


379.5 


vD88 


NO 


2.61 x 10" 10 


0.328 


218.9 


vD88 


S 


6.06 x 10 -10 


0.343 


480.2 


vD88 


H 2 S 


7.40 x 10- 10 


0.344 


594.4 


vD88 



:^uifuiVuXiI(X u i) s~ 



(3) 



where I is the intensity of the radiation field, a is the 
continuum cross section and the line transition is at a wave- 
length X u i from a lower level I to an upper level u with os- 
cillator strength /„; , dissociation efficiency rj u for the upper 
level and occupation fraction xi for the lower level. The lower 
limit of the integral (912 A) is the Lyman cutoff wavelength 
for the unshielded interstellar radiation field. From these 
expressions we have (re-)calculated the unshielded interstel- 
lar photorates usin g cross section and oscillator strength 
data primari ly from van Dishoeck (1988), referring to other 
sources (e.g. iRawlings et~al]|l993l ') for species and reaction 
channels not covered in their database. Where we have one 
set of cross section data for a molecule but multiple dissoci- 
ation channels are list ed in the UMIST Dat abase for Astro- 
chemistry (UMIST06, [Woodall et al]|2007l ). our calculated 
rate is divided between the channels in the same ratio as 
the original UMIST06 rates. We n ote that the h igh resolu- 
tion N2 dissociation line data from ICarterl j 19721 ) has been 
re -binned to a l o wer w avelength resolution in the database 
of Ivan Dishoeck! (|l988T ). As the representative wavelengths, 
we therefore elected to use th e mid-points of the dissocia- 
tion wavelength intervals from lCarterl (fl972). Since no cross 
section data is available for HNC we assign the same values 
as derived for HCN. For the radiatio n field we ado pt the 
standard interstellar radiation field of Drainc ( 1978) , with 
the ex tension to A > 2000 A given bv lvan Dishoeck fc Black! 
(1982). These photorates are presented in Tables [1] and [2] 

The newly calculated direct photorates are in agreement 
with the values in UMIST06 to within a factor of two for all 
but 10 of the 61 reactions for which we had data available. 
Four show rates differing by more than one order of mag- 
nitude; our calculations for the dissociation of Hj" to either 

+ H or H 2 + H+ are each larger by a factor of ~ 400, 
while the ionization of OH is larger by a factor of ~ 12 and 
the dissociation of OH + to + + H is larger by a factor of 
~ 19. 



2.2 Shielding of N 2 and CN by H 2 

The photodissociation of H2 is effected by the absorption 
of Lyman and Werner band photons followed by decay into 
the dissociation continuum (IStecher fc Williams! 119671 ). To 
investigate the possible effects of mutual shielding by these 
lines we note that the range of wavelengths for H2 dissoci- 
ation (844.8 — 1108.5 A) fully o verlaps that for the dissoci- 
ation lines of N 2 (914 - 980 A, Ivan Dishoecklll988l . ICarterl 
1 19721 ) and nearly the whole range for CN photodissociation 
(912 — 1133 A, peaking at 940 A). We are investigating the 
photochemistry well within a PDR (A„ > 2) and the H2 lines 
will be extremely optically thick, to the extent that (a) the 
overlap with individual N2 lines will be extensive and (b) 
there will be significant blanketing of the Lyman/ Werner 
continuum. To include these effects we therefore make the 
extreme assumption that N2 and CN are both fully shielded 
and that the effective rates of photodissociation can be set to 
zero. These assumptions are clearly very simplistic and will 
need addressing using more sophisticated radiative transfer 
techniques in a future study. 



2.3 Shielding by Cl 

To account for shielding by the Cl ionization continuum we 
take a more detailed and accurate approach. First of all, con- 
sidering the situation when the continuum is very optically 
thick, the photoreaction rates are re-calculated using equa- 
tions [2] and [3] but using a lower limit of 1102 A for the in- 
tegral and excluding line contributions for A < 1102 A. This 
implicitly assumes that the radiation field intensity is com- 
pletely attenuated at these wavelengths shortwards of the 
Cl ionization potential; I(A < 1102 A)=0. In Tables Q] and 
[2] we give the resulting shielding factors (Si), which repre- 
sent the ratios of the photoreaction rates without and with 
the contributions from the 912 < A < 1102 A wavelength 
range. Thus, in conditions where the carbon ionization con- 
tinuum is very optically thick, the rate coefficient is just Si 
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Table 2. Photodissociat ion Rates, Carbon Continuum Shielding Factors, Cosmic Ray Induced P hotoreaction Efficiencies an d Cross 
Section Data References; Ivan Dishoeckll 19881 ( http: / /home. strw.leidenuniv.nl/ewine/photo/ 1 vD88). iKulander fc Bottcher|[l978l (KB78) 
and lCarter|[l973 (C72). 



Species 


Photorate / s 


Maximal Carbon 
Shielding Factor, Si 


CR Induced Photoreaction 
Efficiency, pj^f 


Reference 


H+ -» H+ + H 


5.73 x 10- 10 


0.778 


427.7 


vD88 


H+ -> H+ + H 


2.16 x 10~ 12 


0.000 


0.15 


KB78 


H+ -> H 2 + H+ 


2.16 x KT 12 


0.000 


0.15 


KB78 


CO -> O + C 


2.59 x 10 -10 


0.000 


105.0, see text 


vD88 


CO+ -v C+ + O 


1.01 x KT 10 


0.954 


45.7 


vD88 


CH -> C + H 


9.19 x icr 10 


0.975 


465.4 


vD88 


CH+ -> C + H+ 


3.27 x lO" 10 


0.036 


254.7 


vD88 


CH 2 -> CH + H 


5.81 x icr 10 


1.000 


200.4 


vD88 


CH+ -> CH+ + H 


4.53 x 10 -11 


0.461 


66.1 


vD88 


CH+ -> CH + H+ 


4.53 x 1CT 11 


0.461 


66.1 


vD88 


CH+ -> C+ + H 2 


4.53 x 1CT 11 


0.461 


66.1 


vD88 


CH 3 -> CH 2 + H 


1.38 x 10 -10 


1.000 


76.1 


vD88 


CH 3 -7 CH + H 2 


1.38 x icr 10 


1.000 


76.1 


vD88 


CH 4 -> CH 3 + H 


1.89 x 10~ 10 


0.703 


182.8 


vD88 


CH 4 -> CH 2 + H 2 


8.42 x 10~ 10 


0.703 


182.8 


vD88 


CH 4 -> CH + H 2 + H 


1.89 x icr 10 


0.703 


182.8 


vD88 


CH+ -> CH+ + H 2 


2.23 x 10- 10 


0.405 


139.3 


vD88 


CH+ -> CH+ + H 


5.23 x KT 11 


0.405 


139.3 


vD88 


NH -> N + H 


5.02 x 10" 10 


0.794 


219.2 


vD88 


NH 2 — > NH + H 


7.45 x lO" 10 


0.783 


152.8 


vD88 


NH 3 -J- NH + H 2 


4.83 x lO" 10 


0.780 


159.2 


vD88 


NH 3 -> NH 2 + H 


6.82 x lO" 10 


0.891 


148.1 


vD88 


N 2 -» N + N 


2.28 x 10~ 10 


0.000 


184.3 


vD88, C72, see text 


o 2 -> o + o 


7.90 x 10- 10 


0.890 


361.3 


vD88 


o+ -> o+ + o 


3.46 x 10" 11 


1.000 


3.0 


vD88 


OH -> O + H 


3.76 x 10- 10 


0.752 


279.3 


vD88 


OH+ — > 0+ + H 


1.30 x 10" 11 


0.008 


10.7 


vD88 


H 2 — > OH + H 


7.54 x 10" 10 


0.746 


527.0 


vD88 


HCO -> CO + H 


1.11 x lO^ 09 


1.000 


319.8 


vD88 


HCO+ -> CO+ + H 


5.39 x 10~ 12 


0.000 


0.0 


vD88 


TT /~i /~\ , /~( i T T i TT 

H 2 CO — > CO + H + H 


5.79 X 10 


0.844 


132.1 


vD88 


H 2 CO -v CO + H 2 


5.79 x lO" 10 


0.844 


132.1 


vD88 


H 2 CO -+ HCO+ + H + e" 


1.16 x 10" 11 


0.844 


132.1 


vD88 


co 2 -> CO + O 


8.81 x 10~ 10 


0.329 


643.8 


vD88 


CN -> N + C 


2.92 x 10" 10 


0.018 


229.2 


vD88 


HCN -v CN + H 


1.56 x lO" 09 


0.677 


436.9 


vD88 


HNC -v CN + H 


1.56 x 10~ 09 


0.677 


436.9 


vD88, see text 


NO -> O + N 


4.73 x 10~ 10 


0.737 


163.5 


vD88 


HS -> S + H 


9.79 x 10~ 10 


0.947 


271.9 


vD88 


HS+ -v S+ + H 


2.60 x lO" 10 


0.343 


171.9 


vD88 


H 2 S — > HS + H 


1.55 x 10 -09 


0.845 


412.6 


vD88 


H 2 S -> S + H 2 


1.55 x 10~ 09 


0.845 


412.6 


vD88 


cs -> s + c 


9.75 x 10~ 10 


0.973 


1570.0 


vD88 



times the unshielded value. These ratios are calculated for 
all species/channels where cross section data is known. 

We note that the photoionization cross section of atomic 
carbon is oc = 1.6 x 1CP 17 cm 2 and is app roximately in- 
depe ndent of wavelength for A < 1102 A l|van Dishoeckl 
1988) . Cl column densities of up to the order of 10 19 cm -2 
i Tielens &: Hollenbachlll985r ) can be expected in PDRs, lead- 
ing to a wide range of potential optical depths. Taking the 
first approximation that attenuation by dust and molecules 
are separable and the fact that the cross section is uni- 
form with wavelength simplifies the radiative transfer to give 
shielding functions, Oj, for each photoreaction i (analogous 
to those used for CO): 



@i(N c ) = 5' l + (l-S' l )e- <TcJVc , (4) 



where Si is the relevant shielding factor for the optically 
thick limit, given in Tables [1] and [2l and Nc is the carbon 
(Cl) column density. Strictly speaking, the Av-dependence 
(due to continuum absorption by dust; 7 in equation [T]) will 
also depend on the Cl continuum opacity. However, this is a 
complicating second-order effect and we follow the practice 
of previous studies by decoupling the shielding by Cl from 
that by dust absorption, as in equation [T] 
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2.4 Cosmic Ray Induced Photoreactions 

At extinctions greater than a few magnitudes, or if any of the 
various shielding effects described above become important, 
cosmic ray induced photodissociation and photoionization 
processes are significant and must be included in the models. 
The radiation field is generated by seco ndary electrons pro- 
duce d in the cosmic ray ionization of H2 (|Prasad fc Tarafdar] 
19831). The calculation of these rates is as described in 
Sternberg et~ai] (| 19871 ) and lGredel et all (|l989D : 

On the assumption that the total absorption cross sec- 
tion is dominated by dust (rather than the molecular com- 
ponent, which in normal circumstances is true, except for 
H2 and CO) then the photorates are given by: 

R. = T^rPicm-V 1 , (5) 
(1 -u) 

where £ is the cosmic ray ionization rate, n(i) is the 
abundance of species i, and uj is the grain albedo. The cosmic 
ray induced photoreaction efficiency, pi, is given by: 

f <Jj(v)P{v) J 

Pi = ]^^ dV > (6) 

where en is the photoreaction cross section and a g is 
the grain extinction cross section per H-nucleon. The factor 
of 2 in the denominator takes into account the fact that the 
definition of a a is per H-nucleon, whereas £ is defined per 
H 2 molecule (|Woodall et al.ll2007h . We have re-calculated 
the values of pi using the same cross section data as in Sec- 
tion 12.11 and the (high resolution) cosmic ray induced H2 
emission spectrum, P(y) ( Gredel, personal com munication, 
as depicted in Figure 1 of iGredel et al. 1 1 19891 ). This spec- 
trum has been normalised to take account of the various 
transition probabilities and excitations per cosmic ray ion- 
ization. Thus it includes contributions from excitations to 
the B n E+, B' a Z+, C 1 ^, D 1 !^, and D n U u Ry- 

dberg states and to the valence EjF 1 ^ and a 3 Sj" states. 
Also included are excitations into the repulsive 6 3 £i" state 
and the vibrational levels of the ground state (X 1 ^), to- 
gether with cascades to the state. The ratio of the 
H2 populations in J=0 to J=l is taken to be 1 : 3. Results 
of these calculations are also given in Tables [T] and [2] Note 
that we have not re-calculated the cosmic ray induced rate 
for the photodissociation of CO, for which multiple line over- 
lap occurs and the values given in the UMIST06 database 
are reasonably accurate. 

For 24 of the 61 reaction channels, cosmic ray induced 
photoreaction efficiencies are not specified in the UMIST06 
database. Of the remaining 36 (excluding CO), all but 15 
agree to within a factor of 2. Relative to the UMIST06 
database we found the rates for the ionization of NH, the 
dissociation of CN and the dissociation of H2CO to H2 + 
CO are all lower by factors of ~ 12, 25 and 10 respectively. 



3 THE CHEMICAL MODEL 

To investigate the impact of the shielding mechanisms dis- 
cussed above on chemical abundances, we implement them 
in a time- and depth-(A„) dependant chemical model. Since 
we are not attempting to replicate any specific astrophysi- 
cal environment, we adopt physical and chemica l parameters 
similar to the PDR code benchmarking effort of Rol lig et al.l 



Table 3. Models 



Model 


Definition 


X 


n / cm 


Fl 


Standard 


10 


10 3 


F2 


Bright 


10 s 


10 s 


F3 


Dense 


10 


10 5.5 



Table 4. Physical Parameters 



Temperature, T 


50 K 


CR Ionization Rate, f 


5 x 10~ 17 s 


Grain Albedo, lj 


0.5 


X(He) 


10- 1 


X(C) 


KT 4 


X(O) 


3 x 10~ 4 


X(N) 


8 x 10~ 5 


X(Na) 


10~ 6 


X(S) 


10~ 6 



j2007h . These are summarised in Tables [3] and [4] In par- 
ticular, we do not attempt to solve the equations of ther- 
mal balance and the temperature and density (T and n) are 
fixed. We also assume that the PDR is dynamically static 
so that the extinction, A„, is only a function of position. 
The incident radiation field used is the standard interstellar 
radiation field as described above, scaled by the factor \- 

Our chemistry consists of 7 elements in 81 gas phase 
species (Table [5]) linked by 1129 reactions. For the sake of 
simplicity, we do not include any gas-grain chemical re- 
actions. In any case, at the assumed temperature of 50K, 
gas-grain interactions will have negligible significance. The 
elemental helium, carbon and oxygen abundances are the 
same as m iRollig et al l |2007l ) while for nitrogen, sodium 
and sulphur, whic h were absent from t heir models, we take 
abund ances from I Savage et al.l |l992) and lAsplund et al.l 
(2005). Reaction data and rate coefficients are taken from 
the UMIST06 database except where we calculated updated 
photorates as described in SectionfS] We assume a fixed cos- 
mic ray ionization rate of 5.0 x 10 17 s _1 and a grain albedo of 
0.5. The cosmic ray induced photodissociation rate for CO is 
modified by a temper ature dependant factor of (T/300K) 1 ' 17 
l|Woodall et al.l [20071 ) yielding a factor of eight reduction in 
the rate at 50K. For species with no cross section data avail- 
able we use a default cosmic ray induced photoreaction ef- 

Table 5. Chemical species used in the model 



Chemical Species 

H, H^, H , H^ - ' ^3" 
He, He+, Na, Na+, e" 
C, C+, C", CO, CO+, CH, CH+, CH 2 , CH+ 
CH 3 , CH+, CH 4 , CH+, CH+ 
N, N+, NH, NH+, NH 2 , NH+, NH 3 , NH+, NH+, N 2 , N+, N 2 H+ 
O, 0+, O2, O+, OH, OH+, H2O, H2O+, H3O+ 
HCO, HCO+, H 2 CO, H2CO+, CO2, CO+, HCO+ 
CN, CN+, HCN, HCN+, HNC 
NO, NO+, HNO, HNO+, HCNH+, H 2 NC+, HNCO+, H 2 NO+ 
S, S+, HS, HS+, H 2 S, H 2 S+, H 3 S+ 
CS, CS+, C 2 S, C 2 S+, HC 2 S+, HCS, HCS+, H 2 CS+ 
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ficiency of 200; broadly representative of the average value 
for other species. 

The incorporation of the various shielding effects uses 
a number of simplifying approximations to allow for a fo- 
cused study on the effects of H2 and Cl shielding. As in 
other studies we treat the attenuation by dust of the radi- 
ation field as separable from other processes and described 
by the usual exponential dependences on extinction of the 
photorates. Strictly speaking, these dependences should be 
modified to take account of the Cl absorption etc., but such 
variations are relatively minor as compared to the other ef- 
fects reported in this paper. 

F or standard PDR models (e.g. model Fl o fRollig et al.l 
l2007h . the transition from atomic to molecular hydrogen oc- 
curs at A„ ~ 0.1, the Oil— >Cl transition is at A„ ~ 1.0 
and carbon is mostly converted to CO for A„ > 2.0. Deeper 
than this we can assume sufficient columns of H2 and CO 
have formed that they are completely self shielded and the 
effective photodissociation rates for these species are zero. 
For the same reason, when we investigate the effects of the 
inclusion of mutual shielding by molecular hydrogen lines, 
we take the extreme case scenario described in Section [5] 
and set the photodissociation rates of N2 and CN to zero. 
We therefore restrict our investigation to A„ > 2.0 where 
shielding by H2 and CO take this simple form and do not 
attempt to consider the details of the PDR itself, which we 
defer to a later study. Within our model we can therefore 
crudely switch the H2 line mutual shielding on or off and 
assign a fixed carbon column density across the whole ex- 
tinction range so as to investigate the interplay between the 
two mechanisms and their effect on the chemistry. 



4 RESULTS 

Examples of the key results obtained from the model are 
given in Figures [1] and [2] Our model was tested against 
the results for C + , C and CO abundances in model Fl of 
Rolli g et all l|2007h . We are able to reproduce abundances 
of C + , C and CO at positions well within the PDR, while 
observing typical sensitivities to the C/O ratio, the temper- 
ature dependence of the cosmic ray induced CO photodisso- 
ciation rate and the chemical age. For a direct comparison 
we have evolved the chemistry to equilibrium. This typically 
takes ~ 3 x 10 6 years, although we integrate to 100 Myr so 
as to be consistent with their calculations. 

The figures give a clear comparison of the chemical 
structure that is obtained with and without the implemen- 
tation of H2 mutual shielding and carbon continuum shield- 
ing. We present abundances for observable molecules with 
significant abundances and showing abundance variations of 
factors greater than 10 as a result of these two mechanisms. 
The results are given in Figure [T] for the models Fl, F2 
and F3. We also show abundances for C + , C and CO (PDR 
carbon transition species) as well as N2 and CN (species 
shielded by H2) for reference. Probably the most important 
effect on the chemistry is the shielding of N2 , for which the 
photodissociation is completely suppressed by both mecha- 
nisms. Thus, very significant enrichments are visible in all 
three environments and particularly in the range of extinc- 
tions A v ~ 2 — 4. This abundance excess obviously drives 
a more vigorous nitrogen chemistry so that the abundances 



of other nitrogen-bearing species, most notably N2H , NH3 
and CN, are significantly enhanced. The photoionization and 
photodissociation of NH3 are both partially shielded by the 
carbon continuum, while CN is fully shielded by molecular 
hydrogen. As is usually the case for molecular ions, photodis- 
sociation of N2H + is not significant and had been excluded 
from our chemical network. The presence of significant col- 
umn densities of nitrogen-bearing molecules at relatively low 
extinctions is an important result. However, it is also inter- 
esting to note that so few species are significantly enhanced 
by the shielding processes. For example, one might have ex- 
pected to see large increases in HCN and HNC, but while 
their enrichments are tightly correlated they are still limited 
to factors of less than 10. Similarly, despite the strong en- 
richments in NH3, the simpler hydrides NH and NH2 show 
correlated enrichments but of less than an order of magni- 
tude. 

One issue with these particular models is the assump- 
tion that the Cl continuum is very optically thick; with 
t(C) = 10 at A„ = 2. Using a typical relation for A v — 
6.289 x W- 22 N Htotal (|Rollig et all 120071 ) depending on the 
total hydrogen column density Nn total allows us to approx- 
imate the atomic carbon column density. For model Fl, a 
column of roughly 10 16 cm -2 is obtained at A„ = 2, rising 
to ~ 10 17 cm' 2 at A„ = 4, implying a Cl optical depth of 
order unity, but somewhat less than t(C) = 10. The other 
models yield lower values, with the bright model not achiev- 
ing a column of 10 17 cm -2 until A„ ~ 6 and the dense model 
only achieving a maximum column of roughly ~ 10 15 cm -2 . 
However, the A„-N_r/ tota! relationship depends on the gas to 
dust ratio, a quantity which is not only poorly constrained 
but can also vary by orders of magnitude for different astro- 
physical environments. So whilst the results of Figure[T]may 
only be valid beyond a certain extinction for our given A„- 
Njf totai relation, the results become valid at lower A„ as the 
dust to gas ratio falls (since the Cl column density would 
increase for a given A„). We also note that for model F2 
the CoCO transition is pushed much deeper to A„ ~ 8 by 
the more intense radiation field. This probably invalidates 
our assumption of CO being optically thick and perfectly 
self-shielded by A„ = 2. However, one effect of this would 
be to increase the abundance of atomic carbon at intermedi- 
ate depths and hence the optical depth of the Cl continuum 
shielding CO. 

Figure shows the range of enrichments achieved when 
using different levels of shielding in model Fl. The final 
panel shows that as carbon becomes optically thick, the re- 
sulting enrichments are very similar with or without H2 line 
shielding, even though CN is not fully shielded by the car- 
bon ionization continuum. More importantly, the first panel 
shows that in environments where carbon column densities 
are low and only H2 line shielding occurs, significant enrich- 
ments in N2, N2H + and NH3 are still possible. Qualitatively 
similar behaviour is also seen for models F2 and F3. The 
exception is CN, which in model Fl is seen to be coupled 
to the enrichment in atomic carbon through the formation 
channel: 

C + NO -> CN + O. 

Since the H2 photodissociation lines become optically thick 
at low extinctions (A„ ~ 0.1), the effects of H2 mutual 
shielding will always be significant for A„ > 2. We there- 
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Figure 1. Abundances for models with H2 mutual shielding and t(C) = 10 (left), without H2 mutual shielding and t(C) = (centre) 
and relative enrichments between the two (right). Top: Standard (Fl), Middle: Bright (F2), Bottom: Dense (F3). 



fore expect the true chemical enrichment for A v > 2 to be 
some intermediate between these two plots; that is to say, 
the effects of mutual shielding by the H2 photodissociation 
lines will always be present even if the carbon ionization 
continuum is optically thin. 

5 CONCLUSIONS AND DISCUSSION 

We have re-calculated the interstellar photodissociation and 
photoionization rates of a range of interstellar molecules, 



both for direct interstellar and for cosmic ray induced radi- 
ation fields. We have also calculated direct photorates which 
include the effect of partial blocking of the incident radiation 
field by the carbon ionization continuum. 

The specific quantitative results that can be drawn from 
our model are clearly limited by its simplicity. However, us- 
ing approximations for shielding of the incident radiation 
field deep in PDRs by molecular hydrogen and atomic car- 
bon, we present the result of significant (more than an order 
of magnitude) enrichments in the nitrogen-bearing species 
N2, N2H + , NH3 and CN. These results are seen for a range 
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Figure 2. Enrichments from model Fl for a range of H2 and Cl shielding regimes, all relative to the case where there is no H2 mutual 
shielding and t(C) = 0. Note that in the first panel, the curves for N2 and N2H + coincide. 



of physical parameters, including density, radiation field in- 
tensity and carbon column density, but should be most 
prominent in dust poor environments. These species are all 
important observationally in inferring various characteris- 
tics of the PDRs i n which they are detected. For example, 
iBavet et al.l (|2009h have used PDR models to show that 
CN traces the density and relative nitrogen abundance of 
PDR dominated galaxies. First results from the Herschel 
Space Observatory towards W31C have shown nitrogen hy- 
dride (NH3, NH2 and NH) abundances that cannot be ex- 
plained using simple gas-grain or PDR chemical modelling 
ijPersson et al.| [201C0. Other nitrogen-bearing species such as 
HCN and HNC, also believed to trace dense molecular gas, 
show smaller (< 0.5 magnitudes), correlate d enrichments. 
They are often used in bot h galactic (e.g. iRoberts et all 
|201 lT > and extragalactic (e.g. iKamenetzkv et al. fefjlll ) envi- 
ronments to classify the gas chemistry as photon dominated 
or not. It is clear from our work that the shielding by H2 and 
C have a large impact on the abundances of such nitrogen- 
bearing species and both mechanisms need to be included 
in future PDR models. 



apply the shielding factors to more rigorous PDR models 
with self consistently calculated column densities to investi- 
gate the effects on the chemistry in regions of lower visual 
extinction (At, < 2). This is particularly important given 
that our observed enrichments appear more pronounced at 
lower extinctions where molecular shielding by H2 becomes 
more significant relative to dust extinction. It would also be 
useful to do a thorough re-analysis of the cosmic ray induced 
photodissociation rate of CO, which is the dominant source 
of atomic carbon in regions of high extinction (A„ > 4 — 5). 
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Current ongoing studies suggest these mechanisms have 
only a limited effect on column densities in high extinction 
clouds and cannot explain, for example, the high column 
densities of nitr ogen hydrides seen in the interstellar medium 
towards W31C jPersson et al.ll20ldl . In future work we will 
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